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CONFIDENCE  INTERVALS  FOR  THE  RELIABILITY  OF  MULTI -COMPONENT  SYSTEMS 

ABSTRACT 

A  procedure  for  finding  a  confidence  interval  for  the  reliability, 

F,  of  a  multi -component  device  is  presented  which  utilizes  Bernoulli 
test  data  pertaining  to  the  component  parts. 

In  particular,  n^,  Bernoulli  trials  are  carried  out  for  the  1th 
part  of  a  system  built  up  from  k  different  parts.  On  the  basis  of  the 
number  of  observed  failures  X^  (i*l,...,k)  an  interval  estimate  of 
the  reliability  (or  probability  of  functioning)  of  the  system  is  con¬ 
structed. 

In  the  series  case,  with  k  parts,  the  minimum  and  maximum  of  P 

(as  a  function  of  q  's)  is  found  given  the  condition  thac  En  q.  =  I, 

1  tii  *  * 

qi  denoting  the  probability  of  the  i  part  failing  in  a  single  trial. 

The  bounds  are  thus  functions  of  5.  Since  Z  X^  has  expectation  I  and 

is  closely  approximated  by  well-known,  simple  distributions,  a  confidence 

interval  for  t,  and  hence  for  the  bounds  of  P,  and  thus  for  P,  can  be 

found. 

In  the  general  case,  that  is,  for  a  system  having  a  mixture  of 
series  and  parallel'  .onnected  parts,  parametric  minima  and  maxima  for 
P  are  found  using  several  approaches,  and  independent  confidence 
intervals  are  used  to  construct  an  interval  for  P.  A  comparison  of 
the  results  using  this  approach  with  the  results  of  other  procedures 
is  made  for  some  numerical  examples. 
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INTRODUCTION 


The  reliability,  P,  of  a  device  may  be  define!  simply  as  the 
probability  that  it  functions  correctly  under  specified  conditions. 

For  many  devices  an  interval  estimate  of  this  number  P  is  desired.  It 
is  sometimes  impractical,  if  not  impossible,  to  test  many  assembled 
devices,  and  the  suggestion  has  been  made  that  parts  be  tested  sepa¬ 
rately,  and  the  results  combined  to  form  an  interval  estimate.  The 
present  report  indicates  a  means  of  constructing  an  interval  for  P  with 
any  desired  confidence  coefficient,  preferably,  however,  greater  than 
.90.  As  is  the  case  with  all  non-Neyman-shortest  confidence  intervals 
based  on  discrete  variables ,  the  actual  probability  of  coverage  will 
generally  be  somewhat  higher  than  the  selected  coefficient. 

For  any  giver,  device  with  k  distinct  parts,  P  may  be  written  in 
terms  of  q  ,.,,,q  ,  where  q.  denotes  the  probability  of  the  1th  part 
failing.  The  parts  are  assumed  to  function  independently  of  one  another 
within  the  device  in  the  same  way  and  under  the  same  conditions  as 
when  tested  separately.  The  problem  considered  here  is  that  of  finding 
a  confidence  interval  for  P,  given  that  n  Bernoulli  trials  have  been 
conducted  on  the  i  part  (i«l,...,k)  and  that  X^  failures  have  been 
observed  in  the  trials. 

The  approach  of  this  report  is  to  solve  the  problem  first  for  the 
simple  series  case,  and  then  extend  the  results  to  the  general  case. 

For  a  simple  series  device  consisting  of  k  parts,  P  (expressed 
in  terms  of  the  q^)  can  be  minimized  and  maximized  under  the  restriction 
Z  =  I .  Thus  P  will  be  bounded  on  either  side  by  functions  of  i . 

Since  X  =  Z  Xi  has  expectation  i  and  is  approximately  Poisson  (or  bino- 
mially)  distributed,  one  can  easily  find  a  confidence  interval  for  I 
based  on  X,  and  by  appropriate  computations,  for  the  two  bounding  function 
of  i,  and  hence  for  P. 

The  extension  from  the  simple  series  to  the  more  complicated 
devices  follows  from  factoring  P  into  its  series  and  non-series  parts. 
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Both  parts  can  be  parametrically  bounded,  the  series  part  as  outlined 
above  and  the  non-series  part  using  several  approaches.  From  these 
bounds  one  can  write  upper  and  lower  bounds  for  P,  and  by  using  several 
independent  confidence  intervals,  a  confidence  interval  for  P  may  be 
obtained. 

Numerical  examples  are  given  which  illustrate  the  procedure  and 
provide  comparisons  with  the  results  of  some  other  currently  used 
methods . 

The  reader  who  wishes  to  apply  the  results  must  be  willing  to 
perform  some  numerical  computations.  Situations  undoubtedly  exist  for 
which  the  present  approach  is  of  little  use,  and  in  many  other  cases 
the  experimenter' s  supply  of  ingenuity  may  be  heavily  taxed  to  provide 
shortcuts  peculiar  to  the  problem  at  hand.  The  present  discussion  1b 
designed  to  suggest  several  solutions  to  the  problem,  and  to  urge  the 
reader  to  choose  the  best  one  available.  Most  likely  this  choice  cannot 
be  made  without  first  performing  "rial  and  error  calculations. 
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PARAMETRIC  BOUNDS  FOR  P  IN  THE  SIMPLE  SERIES  CASE 

A  very  simple  case  Is  a  device  consisting  of  k  different  parts 

connected  in  series.  Assuming  that  the  parts  function  independently, 

the  reliability  P  equals  JJ(l-q  )  where  q  equals  the  probability 
"fch 

of  the  i  part  failing.  We  can  find  (for  reasons  which  may  not  yet 
be  obvious)  the  range  of  possible  values  that  J|(l-q^)  can  assume 
given  that  £  n^  equals  5,  where  the  are  positive  integers,  and 
I  is  a  positive  fixed  number. 

The  problem  of  finding  the  range  of  P  given  5  can  be  solved  several 

ways,  perhaps  most  easily  by  simple  trial  and  error  calculations.  For 

known  and  fixed  5,  the  minimum  trd  maximum  values  of  P,  expressed  in 

terms  of  functions  of  5  and  the  n^,  may  be  found  by  numerical  trial  and 

error  of  various  values  of  the  q^.  When  P  is  near  one  (say  greater  than 

one  half),  it  is  clear  from  the  algebraic  expansion  of  ITU-  q^ )  that 

1  -  £  q  will  dominate,  and  it  follows  that  the  minimum  and  maximum 

of  P  (at  worst,  approximately)  will  be  attained  when  £  q^,  subject  to 

£  n  q  *  is  respectively  maximized  and  minimized. 

For  example,  if  k=5,  n,-500,  n„*250,  n ,*500,  6*1. 5»  the  condition 

i  d  5 

is  500q^  +  250qg  +  300q^  *  1.5>  and  ^  ^  +  is  minljnized  When 

q1  ■  I.5/3OO,  qg  ■  q,  =  0,  and  maximized  when  q^  •  q^  *  0,  *  I.5/250. 

A  little  arithmetic  shows  that  P  i3  indeed  maximized  and  minimized, 
respectively,  at  these  points.  Similar  results  easily  may  be  seen  to 
nold  for  all  values  of  6  less  than  250,  Hence  for  any  fixed  value  of  6 
in  this  runge, 

1  -  6/250  <  P  <  1  -  6/500 

and  P  may  assume  any  value  in  this  interval,  which  has  width  6/5OO. 

The  general  form  of  the  bounds  for  P  can  be  somewhat  complicated, 
depending  on  the  values  of  the  n^  given.  However,  the  following  two 
general  cases  are  of  greatest  interest:' 

A.  If  n^  =  n  =  £  n^/k  for  i=l, . . . ,k, 

then  1  -  6/H  <  P  <  (1  -  6/nk)k 
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B.  If  0  < 

^  then 

I 


i 


l  n,  -  lunin  n  and  0  <  4  <  min  (rain  n  ,  E  n  -  kmin  n. ) , 

i  i  i  i  1  1  t  1 


1  -  4/mln  n,  <  P  <  1  -  4/max  n . ,  where  i=l,...,k. 
i  1  ”  i  1 
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SAMPLING  AND  USE  OF  THE  BOUNDS 


Suppose  that  the  results  of  n  Bernoulli  trials  are  known  for  the 

1  part  of  a  simple  series  device.  If  X^  denotes  the  number  of 
failures  observed,  is  tinomiaily  distributed  with  parameters  n^  and 
q^.  The  variable  X  =  £  X^  does  not  follow  a  distribution  having  a 
simple  form,  but  it  is  easy  to  verify  that  the  expectation  and  variance 
of  X  ere  I  and  4  -  2  n^q  respectively.  Using  the  method  of 
Lagrange  multipliers,  it  is  also  easily  shown  that  given  4*2  n^q. , 

2  n^q.2  attains  its  minimum  value  when  q^  *  4/2  n^  for  all  i,  and 

hence  the  variance  is  largest  when  all  the  q^'s  are  equal.  But  under 

the  latter  circumstances,  X  is  bincmially  distributed  with  parameters 

2  n  and  4/2  n  .  If  X  were  Poisson  distributed  with  parameter  4, 
i  ‘  1  -if- 

then  the  variance  would  also  be  4.  Writing  the  variance  of  X  under  the 
above  binomial  assumptions  as  V  (X),  and  the  variance  of  X  under 
Poisson  assumptions  as  Vp(X),  the  following  inequality  holds  for  all 
values  of  n^  and  q^,  i=l,...,k: 

V(X)  <  V  (X)  <  V_(X)  or 

4  -  2  n^2  <  |  -  42/2  nt  <  4 . 

Thus  referring  X  to  either  the  Poisson  or  the  appropriate  binomial 
distribution,  one  should  not  be  surprised  to  find  greater  probabilities 
for  the  extreme  tail  values  than  under  the  true  distribution.  If  so, 
when  using  binomial  or  Poisson  confidence  intervals  with  X,  one  may 
expect  them  to  be  conservative  in  the  sense  that  the  probability  of 
containing  the  unknown  parameter  4  will  be  at  least  as  high  as  the 
confidence  coefficient  of  the  binomisl  or  Poisson  intervals.  In  the 
next  section,  comparisons  will  be  presented  concerning  these  proba¬ 
bilities. 

For  the  present,  assuming  that  in  referring  X  to  the  Poisson  or 
binomial  confidence  intervals  one  will  not  be  led  astray,  it  is  a 

*  Were  each  Poisson,  X  would  also  be  Poisson  distributed. 
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simple  matter  tc  look  up  one  or  two-sided  confidence  limits  for  P, 
once  the  appropriate  parametric  bounds  have  been  chosen.  For  example, 
if  the  a-level  two-sided  Poisson  or  binomial  confidence  limits  for  I 
based  on  X  turn  out  to  be  t  ,  t3,  then  whenever  t^  <  £  <  t2,  it  is  also 
true  that  (using  the  case  of  n  =  n  here) 

1  '  V"  <  1  -  S/11  <  P  <  (1  -  l/nk)k  <  (1  -  tj/nk)* 

and  the  resulting  confidence  interval  covers  P  with  probability  at 
least  l-o.  There  will  be  a  growth  in  actual  confidence  coefficient  in 
audition  to  that  due  to  "too-conservative"  intervals  because  P  is  a 
fixed  number  between  two  functions  of  5,  and  the  confidence  interval 
coders  both  of  these  functions. 

In  many  cases,  the  q  will  be  known  to  be  close  to  zero,  and  the 
values  of  n4  will  be  such  that  I  may  be  safely  assumed  to  be  within  the 
range  of  case  B  of  the  previous  section,  or  preferably,  case  A.  There  is 
a  way  of  converting  all  cases  into  the  equal  sample  size  case,  which 
unfortunately  throws  away  part  of  the  available  information  and  subjects 
the  confidence  interval  for  P  to  additional  fluctuation.  This  will  be 
discussed  further  in  the  section  on  complex  systems. 
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THE  BINOMIAL  AND  POISSON  APPROXIMATIONS 

In  this  section  the  distribution  of  X  =  £  is  in  specific 
cases  compared  to  the  distribution  of  X  when  referred  to  the  Poisson 
or  appropriate  binomial  distribution  discussed  above.  A  comparison 
of  the  confidence  coefficients  for  5  under  these  circumstances  is 
also  made. 

For  the  following  cases,  the  density  of  X  has  been  computed  to 
five  decimals: 


(1)  n1  =  n2  =  15,  n^  +  n^  =  5  for  qx  =  0,  .01,  .05,  .1  and 
.2  and  the  Poisson  density  with  parameter  5. 

(2)  n1  =  5,  n2  =  15,  n1<l1  +  n^  =  5  for  =  0,  .05,  .15,  -30 
and  .60. 


TABLE  I 

Density  of  X,  +  Xg  when  X^:  B^,  n-^  =  ng  =  15, 

ni<1l  +  W  =  5  =  *  qi) 

Binomial 
Assumptions 


X 

q^=0  or  .2 

q.j  =  .01  or  ,19 

q1=.05  or  .15 

q^.l 

Poisson  density 

0 

•035lft 

.03646 

,04047 

.04239 

.04979 

1 

.13154 

.13380 

.13908 

.14150 

.14956 

2 

.23090 

.2304b 

.22868 

,227bo 

.22404 

3 

.25014 

.2474 1 

.23952 

.23609 

.22404 

4 

.I87b0 

.18542 

.17949 

.17707 

.16803 

5 

.10318 

.10297 

.10248 

.10230 

.10082 

0 

.04299 

.04389 

.04035 

.04756 

.05041 

7 

.01382 

.01468 

.01705 

.01804 

.02160 

8 

.00346 

.00590 

.00520 

.00576 

.00810 

9 

. 00007 

.OOO65 

.00153 

.00156 

.00270 

10 

.00010 

.00014 

.00029 

.00036 

.00061 

11 

12 

13 

.00001 

.00002 

.00005 

.00001 

.00007 

.00001 

.00022 

.00006 

.00001 

*  This  column  applies  to  Table  II  also,  but  is  not  listed  there. 
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TABIE  II 


Density  of  X^  +  X2  when  X^:  B^,  q^);  n^  =  5,  ng  =  15, 


rlqi  +  n2^2 


=  3 

Binomial 


Assumptions 


X 

V°,  q3S‘2 

o\ 

r— 1 

II 

rov 

O 

II 

< — 1 
a* 

w*15 

*l^-*3,  92=. 1 

V*6,  V0 

0 

.05518 

.05640 

.05876 

.05460 

.01024 

1 

.  13191* 

.13371 

.15680 

.15185 

.07680 

2 

.25090 

.25047 

.22954 

.25200 

.25040 

3 

.25014 

.24753 

.24285 

.25089 

.54560 

4 

.18760 

.0.8554 

.18212 

.18711 

.25920 

5 

.10518 

.10300 

Y-* 

8 

CD 

VJl 

.10240 

.07776 

i/ 

.04299 

.04586 

.04557 

.04272 

7 

.01582 

.01465 

.01601 

.01594 

8 

.00546 

.00588 

.00459 

.00562 

9 

.OOO67 

.00085 

.00108 

.00075 

10 

.00010 

.00014 

.00021 

.00015 

11 

.00001 

.00002 

.00004 

.00002 

In  the  two  tables ,  n^q^  +  n^q^  is  always  three,  and  the  frequency 
functions  may  he  compared  directly  with  the  approximating  binomial  and 
Poisson  frequency  functions.  In  Table  I,  the  maximal  variance  of  X 
occurs  when  q1  =  .1,  and  in  Table  II,  when  q^  *  .15.  It  will  be  noted 
in  both  tables  that  for  the  extreme  tail  values  of  X,  the  approximating 
binomial  and  Poisson  probabilities  are  too  large.  It  is  alsyji  apparent 
that  since  the  probability  that  X  equals  x  is  a  continuous  function  of 
q^,  that  for  all  values  of  q^  (and  hence  q^)  in  both  tables,  the  binomial 
and  Poisson  probabilities  will  be  too  large  in  both  tails,  and  for  more 
values  of  X  in  the  upper  than  in  the  lower  tail.  Also,  in  the  tables, 
the  binomial  and  Poisson  tail  probabilities  are  greater  than  the  corre¬ 
sponding  true  probabilities  for  nearly  the  same  values  of  X  (differing 
at  most  by  one  point),  regardless  of  the  value  of  q^.  The  actual  prob¬ 
abilities  in  both  tails  increase  as  q^  and  q0  approach  equality,  although 
they  are  always  less  than  the  corresponding  Poisson  probabilities, 
suggesting  that  for  arbitrary  k,  the  distribution  of  X  when  all  the 


"  This  is  for  the  true  distribution  of  X  =  X^  +  X^, 
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q^'s  are  equal  is  "closer"  to  the  corresponding  Poisson  distribution 
than  when  twc  q^'s  differ.  This  is  not  surprising,  for  in  case  the  q^'s 
are  all  equal,  X  is  binomially  distributed  with  parameters  £  n^,  £/£  n^, 
and  if  I  ;s  small  relative  to  £  n^,  the  Poisson  approximation  to  the 
binomial  is  known  )o  be  very  good. 

In  this  optimum  case,  wherein  X  is  binomially  distributed,  one  may 
easily  compare  the  tabled  binomial  and  Poisson  densities.  For  a  single 
binomial  variable  X  having  parameters  n,  F  =  l/n,  the  tail  probabilities 
have  been  compared  with  the  corresponding  tail  probabilities  of  the 
Poisson  distribution  (parameter  nF  =  I)  in  the  following  cases: 

(1)  nF=l;  n=5,  F=.2;  n=10,  F=.l;  n=100,  F=.01 

(2)  nF-lOj  n=50,  F=.2;  n=100,  T=.lj  n=10C0,  F*.01 

(5)  nl=50;  n=100,  F=.5i  n-500,  F=.lj  n=1000,  F*.05 

In  all  cases,  exactly  the  same  remarks  as  made  for  the  sum  of  two 
binomial  variables  apply,  with  the  obvious  additional  observation  that 
as  n  increases,  the  tail  probabilities  of  the  binomial  distribution 
approach  the  corresponding  Poisson  probabilities.  Both  the  binomial 
and  Poisson  distributions  tre  asymptotically  normally  distributed 
(that  is,  as  n  and  i,  respectively,  increase)  and  hence  if  X  is  large, 
one  may  expect,  the  difference  between  the  two  sets  of  confidence  limits 
to  be  relatively  small. 

In  any  particular  problem,  the  preceding  indicates  that  one  will 

achieve  narrower  confidence  bounds  for  i  and  hence  for  P,  by  using  the 

binomial  rather  than  Poisson  limits.  When  i  is  small  and  £  n,  is 

i 

relatively  large,  the  difference  between  the  two  approximating  dis¬ 
tributions  becomes  quite  small.  When  £  is  large  and  X  turns  out 
to  be  small,  tne  binomial  confidence  limits  are  easily  and  quite 
accurately  approximated  by  the  Poisson  limits.  In  the  following  tables 
and  remarks,  comparisons  aie  nade  of  the  Poisson  and  binomial  intervals 
for  certain  small  values  of  n  =  £  n^,  a,  and  X  <  }0,  For  large 
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values  of  n,  the  Poisson  limits  will  be  much  closer  to  the  binomial 
limits  than  for  those  values  tabled. 

These  tables  have  been  drawn  from  three  widely  available  sources: 

Tables  III  and  VI  are  from  the  binomial  graphs  and  Poisson 
confidence  limit  tables  in  3iometrlka  Tables  for  Statisticians , 

Volume  I,  edited  by  E.  S.  Pearson  and  H.  0.  Hartley.  They  are  based  on 
the  equal-tails  approach,  which  chooses,  for  a  given  value  of  X  (say  c) 
the  values  of  a  Poisson  parameter  such  that,  in  the  two-sided  case, 

Z  e"-  xVi’.  =  n/2  =  Z  e"*  X1/!'. 
i=c  1=0  •• 

where  (\,  4)  is  to  be  the  a-level  confidence  interval.  The  two-sided 
intervals  are  tabled  for  Poisson  X  =  0(1)30(5)50,  a  =  .10,  .05,  .02, 

.01,  .00?,  and  binomial  X  <  n,  for  n  <  1000,  a  =  .05,  .01. 

Tables  IV  and  VII  are  based  on  tatles  published  by  Crow  and  Gardner 
in  Biometrika,  as  follows: 

(1)  Crow,  E.  L.  Confidence  Intervals  for  a  Proportion,  Biometrika, 
Volume  43  (1956),  pages  423-435.  (For  binomial  X  <  n, 

n  =  1(1)30,  a  =  .10,  .05,  .01) 

(2)  Crow,  E.  L.  and  Gardner,  R.  S.  Confidence  Intervals  for  the 
Expectation  of  a  Poisson  Variable,  Biometrika,  Volume  46 
(1959),  pages  441-453.  (For  X  =  0(1)300,  a  =  .20,  .10,  .05, 
.01,  .001)  The  system  used,  described  in  detail  in  (l),  is 
optimum  in  a  geometrical  sense  and  generally  yields  bounds  of 
less  width  than  does  the  equal-tails  method. 

Tables  V  and  VIII  are  based  on  tables  published  by  Blyth  and 
Hutchinson  in  Biometrika,  as  follows: 

(1)  Tables  of  Neyraan-Shortost  Unbiased  Confidence  Intervals  for 
the  Binomial  Parameter,  Biometrika,  Volume  47  (i960),  pages 
381-391.  (For  X  <  n,  n  =  2(1)24(2)50,  a  =  .05,  .01) 

(2)  Tables  of  Neyman-Shortest  Unbiased  Confidence  Intervals 
for  the  Poisson  Parameter,  Biometrika,  Volume  48  (1961), 


pages  191-19*1.  (For  X  =  0(l)?90,  a  =  .09,  .01) 

The  optimum  property  of  these  intervals  may  he  described  as 
follows:  Among  ail  unbiased  a-level  confidence  intervals,  the  tabled 
intervals  uniformly  minimize  the  probability  of  covering  false  values. 
An  unbiased  interval  A  is  defined  such  that  if  P^  denotes  probability 
when  the  distribution  parameter  is  0,  that  P  (GeA)  <  1  -  a  foi  all  G, 

V  ^ 

and  P  (©'  eA)  <  P_(GeA)  for  all  G,  G* . 

u  — •  y 

In  the  binomial  portions  of  each  table,  the  entries  are  nc^,  nc^, 
where  c1 ,  cg  are  the  tabled  a-level,  two-sided  confidence  limits  for  p. 

When  using  the  Neyman-shortest  tables,  it  is  necessary  to  choose 
■i  random  number  between  zero  and  one  and  add  it  to  X  and  then  read  the 
appropriate  entry. 
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TABLE  III 

* 

Two-sided  Poisson  vnd  Binomial  Confidence  Limits: 
Equal  Tails,  a  =  .05 


X 

n 

*10 

n= 

.20 

n= 

30 

Poisson 

0 

.00 

3.1 

.00 

3.4 

.00 

3.5 

.000 

3.69 

1 

.02 

4.4 

.02 

5.0 

.06 

5.2 

.0253 

5.57 

2 

.25 

5.6 

.25 

6.3 

.22 

6.6 

.242 

7.22 

5 

.68 

6.5 

.65 

7.5 

.63 

8.0 

.619 

8.77 

4 

1.2 

7-4 

1.2 

8.7 

1.1 

9.3 

1.09 

10.24 

5 

1.9 

8.1 

1.7 

9.8 

1.7 

10.4 

1.62 

II.67 

u 

2.6 

8.8 

2.4 

10.8 

2.3 

11.6 

2.20 

13.06 

J 

3.5 

9.3 

3.0 

11.8 

3.0 

12.  T 

2.81 

14.42 

8 

4.4 

9.8 

3.8 

12.8 

3.6 

13.7 

3.45 

15.76 

9 

5.6' 

10.0 

>.6 

13.7 

4.4 

14.8 

4.12 

17.08 

10 

6.9 

10.0 

5.4 

14.6 

5.2 

15.8 

4.80 

18.39 

11 

6.3 

15.4 

5.9 

16.8 

5.49 

19.68 

12 

/  .2 

16.2 

6.8 

17.8 

6.20 

20.96 

13 

8.2 

17.0 

7.6 

18.8 

6.92 

22.23 

14 

9.2 

17.6 

8.4 

19.6 

l  .63 

25.49 

15 

10.2 

18.3 

9.4 

20.6 

8.40 

24.74 

l6 

11.3 

18.8 

10.4 

21.6 

9.15 

pc  <>A 
*-✓ . 

if 

12.5 

19.4 

11.2 

22.4 

9.90 

27.22 

18 

13.7 

19.8 

12.2 

2,’.2 

IO.67 

28.45 

1? 

15.0 

20.0 

13.2 

24.1 

11.44 

29.67 

20 

16.6 

20.0 

14.2 

24 .8 

12.22 

30.89 

21 

15.2 

25.6 

13.00 

32.10 

22 

16.3 

26.4 

13.79 

33.31 

23 

17.3 

27.0 

14.58 

34.51 

24 

18.4 

27. 1 

15.38 

35.71 

25 

19.6 

28.3 

16.18 

36.90 

2o 

20.7 

28.9 

16.98 

38.10 

2 1 

22.0 

29.4 

17.79 

39.28 

28 

23.4 

29.8 

l8.6l 

40.47 

29 

24.8 

29.9 

19.42 

41.65 

30 

26.5 

30.0 

20.24 

42.83 

*  The 

above 

binomial 

confidence  limits 

are  for 

np. 

TABLE  IV 


* 


Two-Sided  Limits: 

Crow  and  Gardner,  a  =  . 

>05 

X 

n=! 

5 

n=10 

n= 

>20 

n= 

=30 

Poisson 

0 

.000 

2.5 

.000 

2.67 

.000 

2.86 

.000 

3.00 

.0 

3.285 

1 

.050 

3.28 

.05* 

3.97 

.06 

4.44 

.06 

4.89 

.051 

5.323 

2 

.380 

4.06 

.37 

6.03 

.36 

5.86 

.36 

6.15 

.355 

6.686 

3 

•  9^5 

4.62 

.87 

6.19 

.84 

7.02 

.84 

7.32 

.818 

8.102 

4 

1.72 

4.95 

1.50 

7.33 

1.42 

8.22 

1.41 

8.76 

1.366 

9.598 

5 

2.50 

5.00 

2.22 

7.78 

2.08 

9.34 

2.04 

9.72 

1.970 

11.177 

6 

2.67 

b.50 

2. Go 

10.7 

2.73 

10.9 

2.613 

12.817 

7 

3.81 

9.13 

2.86 

11.8 

3.00 

12.1 

3.285 

13.765 

8 

3.97 

9.63 

4.18 

13.0 

3.93 

13.2 

3.285 

14.921 

9 

6.03 

9.95 

4.44 

14.1 

4.89 

14.3 

4.460 

16.768 

10 

7.33 

10.0 

5.86 

l4.l 

5.25 

15.7 

5.323 

17.633 

11 

5.86 

15.6 

6.15 

16.8 

5.323 

19.050 

12 

7.02 

15.3 

7.08 

17.9 

6.686 

20.335 

13 

8.22 

17.1 

7.32 

19.1 

6.686 

21.364 

14 

9.34 

17.2 

8.76 

20.3 

8.102 

22.945 

15 

10.7 

17.9 

9.72 

20.3 

8.102 

23.762 

lb 

11.8 

18.6 

9.72 

21.2 

9.598 

25.400 

17 

13.0 

19.2 

10.9 

22.7 

9.598 

26.306 

l8 

14.1 

19.6 

12.1 

22.9 

11.177 

27.735 

19 

15.6 

19.9 

13.2 

23.8 

lJ-.l  (< 

28.966 

20 

17.1 

20.0 

14.3 

24.8 

12.817 

30.017 

21 

15.7 

25.1 

12.017 

31.675 

22 

16.8 

26.1 

13.765 

32.277 

23 

17.9 

27.0 

14.921 

34.048 

24 

19.1 

27.3 

14.921 

34.665 

25 

20.3 

28.0 

16.768 

36.030 

26 

21.2 

28.6 

16.77 

37.67 

27 

22.7 

29.2 

17.63 

38.16 

28 

23.8 

29.6 

19.05 

39.76 

29 

25.1 

29.9 

19.05 

40.94 

30 

27.0 

30.0 

20.33 

41.75 

4~~The  above  ’binomial  confidence  limits  are  for  np. 


TABLE  V 

* 

Two-Sided  Neyman-Shortest,  Limits:  a  =■  .0? 


X+Y 

n= 

=5 

n= 

=10 

n=20 

n=30 

Poisson 

,0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.00 

.00 

.  ’ 

.0 

2.2 

.0 

2.5 

.0 

2.6 

.0 

2.7 

.0 

2.8 

1.0 

.0 

2.6 

.0 

3.0 

.0 

3.2 

.0 

3.3 

.0 

3.5 

1.5 

.0 

5.4 

.0 

4.0 

.0 

4.3 

.0 

4.4 

.0 

4.7 

-'.0 

.0 

3.8 

.1 

4.4 

.0 

4.8 

.0 

5.1 

.0 

5.4 

2.5 

p 

4.3 

.1 

5.2 

.1 

5.7 

.0 

5.8 

.1 

6.4 

5.0 

.5 

4.5 

.4 

5-6 

.4 

6.2 

.3 

6,6 

.4 

7.1 

0.5 

•7 

4.8 

.6 

6.3 

•  5 

7.1 

.6 

7.4 

.5 

8.0 

'+.0 

1.2 

5.0 

1.0 

6.6 

.8 

7.6 

•  9 

7.8 

.8 

8.6 

4.5 

l.b 

5.0 

1.2 

7.2 

1.1 

8.3 

1.0 

8.6 

1.0 

9.5 

5.0 

2.4 

5.0 

1.7 

7.5 

1.6 

8.8 

1.5 

9*3 

1.4 

10.1 

S5 

2.0 

8.0 

1.7 

9*3 

1.6 

9-9 

1.6 

10.9 

0 .0 

"'5 

8.3 

2.2 

9.8 

2.1 

10.5 

2.0 

11.5 

6.5 

-.8 

8.8 

2.4 

10.4 

2.4 

11.1 

2.2 

12.3 

7.0 

5.4 

9.0 

5.0 

10.8 

2.7 

11.4 

2.6 

12.9 

7.5 

5.7 

9.4 

3.2 

11.4 

3-0 

12.3 

2.9 

13.7 

P.0 

4.4 

9.6 

3.6 

11.8 

3.6 

12.6 

3.3 

14.3 

8,s 

4.8 

9.9 

4.0 

12.4 

3.9 

13.2 

3.6 

15.0 

;..b 

5.6 

9.9 

4.6 

12.8 

4.2 

13.8 

5.9 

15.6 

r'.5 

6.0 

10.0 

4.8 

13.4 

4.5 

14.4 

4.2 

16.3 

10.0 

7.0 

10.0 

5.4 

15.8 

5.1 

14.7 

4.6 

16.9 

*  The  above  "binomial  confidence  limits  are  for  np. 
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TABLE  V 


Two-Sided  Neyman-Shortest  Limits:  a  =  .05 


(Continued) 

X+Y 

n=5  n=10  n=20 

n=30 

Poisson 

11 

6.2 

14.6 

6.0 

15.9 

5.3 

18.2 

12 

7.2 

15.4 

6,6 

16.8 

6.0 

19.5 

13 

8.2 

16.4 

7.5 

18.0 

6.6 

20.8 

14 

9.2 

17.0 

8.4 

18.9 

7.5 

22.1 

15 

10.2 

17.8 

9.3 

19.8 

8.2 

23.3 

16 

11.2 

18.4 

10.2 

20.7 

9.0 

24,6 

17 

12.4 

19.2 

11. 1 

21.6 

9.7 

25.8 

18 

13.8 

19.6 

12,0 

22.5 

10.5 

27.1 

19 

15.2 

20.0 

13.2 

25.4 

11.3 

28.5 

20 

16.8 

20.0 

14.1 

24.0 

12.0 

29.5 

21 

15.3 

24.9 

12.8 

30.7 

22 

lb. 2 

25.8 

13.6 

31.9 

23 

17.4 

26.4 

14.4 

33.1 

24 

18.6 

27.3 

15.2 

34.3 

25 

19.5 

27.9 

16.0 

35.5 

26 

• 

20.7 

28.5 

16.8 

36.7 

27 

22.2 

29.1 

17.6 

37.9 

28 

25.4 

29.7 

18.4 

39.1 

29 

24.9 

30.0 

19.3 

40.3 

30 

26.7 

30.0 

20.1 

41.5 

The  above  binomial  confidence  limits  are  for  np. 
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TABLE  VI 

# 

Two-Sided  Equal -Tails  Limits:  a  =  .01 


X 

n= 

=10 

n= 

20 

n=30 

Poisson 

0 

.00 

4.1 

.00 

4.7 

.00 

4.9 

.000 

5.30 

1 

.00 

5.4 

.00 

6.2 

.00 

6,6 

.00501 

7.43 

2 

.12 

6.5 

.12 

7.8 

.08 

8.2 

.10.3 

9.27 

3 

•  37 

7.4 

.40 

8.9 

.38 

9.6 

.338 

10. '98 

4 

.77 

8.1 

•  74 

10.1 

.75 

10.9 

.672 

12.59 

5 

1.3 

3.7 

1.1 

11.2 

1.1 

12.2 

1.08 

14.15 

6 

1.9 

9-2 

1.7 

12.2 

1.6 

13.3 

1.54 

15.66 

7 

2.6 

9.6 

2.3 

13.1 

2.6 

14.4 

2.04 

17.13 

8 

3o 

9.9 

2.9 

14.0 

2.3 

15.4 

2.57 

18.58 

9 

4.6 

10.0 

3.6 

14.8 

3.4 

16.5 

3.13 

20.00 

10 

5.9 

10.0 

4.4 

15.7 

4.1 

17.5 

3.72 

21.40 

11 

5.2 

16.4 

4.7 

18.4 

4.32 

22.78 

12 

6.0 

17*1 

5.6 

19.4 

4.94 

24.14 

13 

6.9 

17.7 

6.3 

20.2 

5.58 

25.50 

14 

7.3 

18.3 

7.1 

21.2 

6.23 

26.84 

15 

8.8 

18.9 

8,0 

22.1 

6.89 

28.16 

16 

9.9 

19.3 

8.8 

22.9 

7.57 

29.48 

17 

11.1 

19.6 

9.8 

23.7 

8.25 

30.79 

l8 

]?.? 

19.9 

10.6 

24.4 

8.94 

32.09 

19 

13.3 

20.0 

11.6 

25.3 

9.64 

35.38 

20 

15.3 

20.0 

12.5 

25.9 

10.35 

34.67 

21 

13.5 

26.6 

11.07 

35.95 

r.'2 

14.6 

27.2 

11.79 

37.22 

'  1  V 

i  c  r, 

•  - 

?7.4 

12.52 

38.48 

24 

16.7 

28.4 

13.25 

39.74 

?5 

17.8 

28.9 

14.00 

41. OC 

26 

19.1 

29.2 

14.74 

42.25 

27 

20.4 

29.6 

15.49 

43.50 

28 

21.8 

29.9 

16.24 

44.74 

^9 

23.4 

30.0 

17.00 

45.98 

30 

25.1 

30.0 

17.77 

47.21 

*  The 

above 

binomial 

confidence  limits 

are  for 

np. 

22 


TABLE  VII 

* 

Two-Sided  Limits:  Crow  and  Gardner,  a  *  .01 


X  n=5  n=10  n=20  n=30  Poisson 


0 

.000 

3.01 

.00 

3.12 

.00 

4.18 

.00 

4.53 

.000 

4.771 

1 

.010 

3.89 

.01 

5.12 

.02 

5.86 

.00 

6.18 

.010 

6.914 

2 

.165 

4.47 

.16 

6.24 

.16 

7.50 

.15 

7.68 

.149 

8.727 

3 

•  550 

4.84 

.48 

7.03 

.46 

8.48 

.45" 

9.30 

.436 

10.473 

4 

1.11 

4.99 

.93 

7.82 

.88 

10.0 

.84 

10.4 

.823 

12.347 

5 

1.99 

5.00 

1.50 

G.50 

1.38 

11.5 

1.35* 

11.6 

1.279 

13.793 

6 

2.18 

9.07 

1.96 

12.0 

1.89 

12.9 

I.785 

15.277 

7 

2.97 

9.52 

2.58 

12.7 

2.49 

14.1 

2.330 

16.001 

8 

3.76 

9.84 

3.26 

14.1 

3.12 

15.2 

2.906 

16.362 

9 

4.88 

9.99 

4.00 

14.5 

3.81 

16.1 

3.507 

19.462 

10 

6.24 

10.0 

4.18 

15.8 

4.53 

17.1 

4.130 

20.676 

11 

5.48 

16.0 

4.53 

18.4 

4.771 

22.042 

12 

5.86 

16.7 

5.94 

19.6 

4.771 

23.765 

13 

7.26 

17.4 

6.18 

20.1 

5.829 

24.925 

14 

7.98 

18.0 

7.47 

20.8 

6.668 

25.992 

15 

8.48 

18.6 

7  »C8 

22.3 

6.914 

27.718 

16 

10.0 

19.1 

9.24 

22.5 

7.756 

28.852 

1/ 

11.5. 

19.5 

9.87 

23.8 

8.727 

29.900 

18 

12.5 

19.8 

10.4 

24.1 

8.727 

31.839 

19 

14.1 

20.0 

11.6 

25.5 

10.009 

32.547 

20 

15.8 

20.0 

12.9 

25.5 

10.473 

34.183 

21 

13.9 

26.2 

11.242 

35.204 

22 

14.8 

26.9 

12.347 

36 . 544 

23 

15.9 

27.5 

12.347 

37.819 

24 

17.1 

28.1 

13.793 

38.939 

25 

18.4 

28.6 

13.793 

40.373 

26 

19.6 

29.2 

15.28 

41.39 

27 

20.7 

29.6 

15.28 

42.85 

28 

22.3 

29.8 

16.80 

45.91 

29 

23.8 

30.0 

16.80 

45.26 

30 

25.5 

30.0 

18.36 

46.50 

*  The  above  binomial  confidence  limits  are  for  np. 
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TABLE  VIII 


* 

Two-Sided  Neyman -Shortest  Limits:  a  =  ,01 


X+Y 

n= 

■b 

n= 

>10 

n=20 

n»30 

Poisson 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

.0 

•  5 

.0 

3.1 

.0 

3.6 

.0 

4.0 

.0 

4.2 

.0 

4.4 

1.0 

.0 

3.4 

.0 

4.1 

.0 

4.6 

.0 

4.8 

.0 

5.1 

1.5 

.0 

4.1 

.0 

5.1 

.0 

5.7 

.0 

6.2 

.0 

6.6 

2.0 

.0 

4.3 

.0 

5*5 

.0 

6.2 

.0 

6,6 

.0 

7.2 

2.5 

.0 

4.7 

.0 

6.2 

.0 

7.3 

.0 

7.6 

.0 

8.4 

5.0 

.2 

4.8 

.2 

6.6 

.2 

7.8 

.3 

8.1 

.1 

9.1 

3-5 

.5 

5.0 

.2 

7.2 

.2 

8.5 

.3 

9.2 

.2 

10.2 

4.0 

*7 

5.0 

.5 

7.5 

.4 

9.0 

.6 

9.6 

.4 

10.8 

4.5 

.9 

5.0 

.7 

8.1 

.6 

9-7 

.6 

10.4 

.5 

11.8 

5.0 

1.6 

5.0 

1.1 

8.3 

1.0 

10.2 

.9 

10.8 

.8 

12.4 

5.5 

1.2 

8.8 

1.1 

10.9 

1.0 

11.6 

1.0 

13.4 

6.0 

1.7 

8,9 

1.4 

11.2 

1.5 

12.0 

1.3 

14.0 

6.5 

1.9 

9.3 

1.6 

11.8 

1.5 

12.9 

1.4 

14.9 

7.0 

2.5 

9.5 

2.0 

12.2 

2.1 

13.2 

1.8 

15-5 

7.5 

2.8 

9-8 

2.2 

12.8 

2.1 

14.1 

2.0 

16.4 

8.0 

3.4 

9.8 

2.8 

15.2 

2.7 

14.4 

2.3 

17.0 

8.5 

3.3 

10.0 

3.0 

13.8 

2.7 

15.0 

2.5 

17.8 

9.0 

4.5 

10.0 

5.4 

14.0 

3.3 

15.6 

2.9 

18.4 

9.5 

4.9 

10.0 

3.8 

14.6 

3.6 

16.2 

3.1 

19.2 

10.0 

5.9 

io.o 

4.2 

15.0 

3.9 

16.5 

3.5 

19.8 

*  The  above  binomial  confidence  limits  are  for  np. 
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TABLE  VIII 

Two-Sided  Neyman-Shortest  Limits:  a  »  .01 
(Continued) 


X+Y  n-5  n«10 


n=20 


n»30  Poisson 


11 

5.0 

15.8 

4.8 

17.4 

4.1 

21.2 

12 

6.0 

16.6 

5.4 

18.6 

4.8 

22.6 

15 

6.8 

17.2 

6.3 

19.5 

5.4 

24.0 

14 

7.6 

l8.o 

6.9 

20.4 

6.0 

25.3 

15 

8.8 

18.6 

7.8 

21.3 

6.7 

26.7 

16 

9.8 

19.0 

8.7 

22.2 

7.4 

28.0 

17 

11.0 

19.6 

9.6 

23.1 

8.1 

29.5 

18 

12.2 

19.8 

10.5 

23.7 

8.8 

30.6 

19 

13.8 

20.0 

11.4 

24.6 

9-5 

31.9 

20 

15.4 

20.0 

12.6 

25.2 

10.2 

33.2 

21 

15.5 

26.1 

10.9 

34.5 

22 

14.4 

26.7 

11.6 

35.8 

23 

15.6 

27.3 

12.4 

37.0 

24 

16.8 

27.9 

13.1 

38.3 

25 

18.0 

28.5 

13.8 

39.6 

26 

19.2 

29.1 

14.6 

40.8 

27 

2U.4 

29.4 

15.3 

42.1 

28 

21.9 

29.7 

16.1 

43.3 

29 

23. 4 

30.0 

16.8 

44.6 

30 

25.2 

30.0 

17.6 

45.8 
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It  is  apparent  from  visually  comparing  the  Poisson  limits  with 
the  corresponding  binomial  limits,  that  except  for  a  few  cases  (a  total 
of  16),  when  a  =  .0?  nr  .01,  the  former  completely  cover  the  binomial 
limits  for  n  =  10,  20,  30  and  in  four  of  the  tables,  for  n  ~  5*  In 
Tables  TIT  and  VI  the  only  discrepancies  occur'  for  X  <  2,  vbp*-e  the 
graphed  values,  multiplied  by  n,  may  be  necessarily  too  imprecise  to 
compare  meaningfully  with  the  Poisson  readings.  In  Tables  IV  and  VIII, 
for  the  discrepancies  occurring  when  X  <  2,  again  impreciseness  of  the 
binomial  tables  may  be  responsible.  For  the  larger  values  of  X  the 
behavior  of  both  intervals  is  somewhat  erratic,  and  this  plus  the  fact 
that  both  the  tabled  Poisson  and  binomial  intervals  generally  cover 
the  unknown  parameters  with  probability  greater  than  .93  or  .99# 

Indicates  that  an  occasional  slight  non-inclusion  for  a  particular 
value  cf  X  is  not  serious.  ?ue  procedure  of  referring  to  the  binomial 
or  Poisson  limits  will  still  be  expected  to  yield  confidence  coefficient 
very  near  to,  if  not  greater  than  1  -  a. 

In  all  the  tables,  as  n  increases,  the  end-points  uniformly  expand 
towards  the  corresponding  values  in  the  Poisson  case.  For  a  =  .05  or 
.01,  the  fit  generally  seems  to  be  closer  for  the  Neymnn-shortest 
intervals,  which  is  not  surprising  since  it  is  the  only  one  of  the 
three  methods  that  exactly  attains  its  confidence  coefficient  for  all 
parameter  values. 

For  values  of  X  and  n  other  tlmm  those  tabled,  several  approximations 
to  the  binomial  and  Poisson  distributions  exist,  a  few  of  which  appear 
in  Biometrika  Tables  for  Statisticians. 

As  a  check  on  the  inclusion  properties  of  the  binomial  and  Poisson 
approximations  one  may  for  specific  examples,  involving  known  parameters, 
compute  the  true  probabilities  of  covering  P  using  the  appropriate 
binomial  or  Poisson  confidence  limits  (a-level).  If  the  Poisson  limits 
include  the  binomial  limits,  which  in  turn  are  "too  wide"  for  P,  then 
one  would  expect  in  most  cases  that  the  Poisson  intervals  would  have 
the  highest  confidence  coefficient,  followed  by  the  "binomial" 
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confidence  coefficient,  both  of  which  would  be  expected  to  be  greater 
than  1  -  a.  Calculations,  based  on  the  distribution  of  X  under  the 
various  values  of  appearing  in  Tables  I  and  II,  were  made  to 
obtain  the  probabilities  that  the  Poisson  limits  and  the  appropriate 
binomial  limits  (a-level)  cover  +  n^.  These  probabilities 
are  given  in  Tables  IX  and  X. 
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TABLE  n 


Probability  of  the  Poisson  Limits  Covering  n^  +  n^  =  5 


nl  =  n2 

=  15  a  =  . 

.05 

Interval 

*1  *  0 

^  =  .01 

qx  =  .05 

qx  =  .10 

Poisson 

Equal-Tails 

•99575 

.99509 

.99512 

.99221 

.98809 

Crow  and  Gardner 

.98195 

.98041 

•  97*507 

.57417 

.96649 

Neyman-Shortest 

•  970 

.969 

.964 

.961 

.952  (-95) 

nl  =  n2 

=  15  a  =  , 

.01 

Interval 

<h.m0 

9-^  =  .01 

»  .05 

<ix  =  .10 

Poisson 

Equal-Tails 

.99921 

.99899 

.99852 

•  99757 

.99619 

Crow  and  Gardner 

.99921 

.99899 

.99852 

.99797 

.99619 

Neyman-Shortest 

•  994 

.994 

.995 

.992 

.990  (.99) 

ni  = 

n2  =  15  a 

■  .05 

Interval 

Gj  =  0 

K'N 

O 

• 

II 

qx  =  .15 

qx  -  .50 

q^  0  .60 

Equal-Tails 

•  99575 

.99514 

.99408 

.99549 

1.00000 

Crow  and  Gardner 

.98195 

.98051 

.97807 

.98155 

1.00000 

Neyman-Shortest 

•  970 

.969 

.966 

.971 

.994 

nl  “  5’ 

n2  =  15  a 

=  .01 

Interval 

ql  =  0 

11  =  *05 

q.x  =  .15 

qx  =  .50 

q-j  *  .60 

Equal -Tails 

.99921 

.99902 

.99867 

.999H 

1.00000 

Crow  and  Gardner 

.99921 

.99902 

.00867 

.99911 

1.00000 

Neyman -Shortest 

.994 

.994 

•  99  4 

•  994 

.999 
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TABLE  X 


Probability  of  the  Binomial  Limits  Covering  a  3 


nl  =  n2 

=  15  a  = 

.05 

Interval 

»!-0 

qx  =  .01 

qx  =  .05 

qL  =  .10 

Equal-Tails 

.99575 

.99509 

.99312 

.99221 

Crov  and  Gardner 

.99575 

.99509 

.99312 

.99221 

Neyman-Shortest 

.962 

.960 

.954 

.951  (.95) 

nl  =  n2 

=  15  a  = 

.01 

Interval 

*1  =  ° 

-  -01 

q^  =  .05 

qx  -  .10 

Equal-Tails 

•99921 

.99899 

.99832 

.99797 

Crow  and  Gardner 

.99575 

•99509 

.99312 

.99221 

Neyman-Shortest 

•995 

•  995 

.991 

.990  (.99) 

nx  “  5, 

ng  -=15  a 

-  .05 

Interval 

»l“° 

ql  “  *0-! 

q-L  =  .15 

qx  --  .30 

q^  =  .60 

Equal-Tails 

.99575 

•99514 

.99406 

.99549 

1.00000 

Crow  and  Gardner 

.96037 

.93071+ 

•95532 

.96089 

.98976 

Neyman-Shortest 

.956 

•  953 

•949( .95) 

.956 

•  992 

nl  = 

ng  =  15  a 

=  .01 

luterv&i 

o,  -  0 

a- 

qx  =  .03 

qx  =  .15 

qx  *  .30 

q^  =  ,60 

Equal-Tails 

•99921 

.99902 

.99867 

.99911 

1.00000 

Crov  and  Gardner 

•99575 

.99514 

.99408 

.99549 

1.00000 

Neyman-Shortest 

-991 

.991 

•990(.99) 

•991 

.998 
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Again,  one  ir&y  notice  that  the  smallest,  probabilities  occur  when 
=  .1  and  q1  =  =  .15,  where  the  variance  of  X  is  largest. 

Also,  the  Neyman-shortest  intervals  have  the  smallest  probabilities 
of  all  three  systems,  in  most  cases,  which  is  due  mainly  to  the  exact¬ 
ness  of  these  intervals.  Owing  to  the  number  of  significant  decimals 
published  in  tables  of  the  latter  intervals,  the  probabilities  of 
coverage  may  be  in  error  by  a  few  places  in  the  last  decimal.  For 
example,  in  the  Poisson  case,  for  a  ■  .05,  =  np  “  ^5/  the  computed 

probability  using  linear  interpolation  turns  out  to  be  .952  whereas 
the  true  probability  is  exactly  .95.  The  corresponding  probability 
when  a  =  .01  is  correct  to  three  decimals. 

When  confronted  with  a  value  of  X  resulting  from  lengthy 
experimental  testing,  the  procedure  adopted  might  be  to  select  from 
the  available  tables  the  limits  which  are  narrowest  for  the  particular 
value  of  X  observed,  and  not  to  choose  the  bounds  with  respect  to 
any  optimal  mathematical  property.  Such  a  procedure  defines  a  kind 
of  composite  interval  about  which  little  may  confidently  be  said 
concerning  the  probability  of  coverage.  Preferably,  a  aingls  type  of 
interval  will  be  decided  upon  in  advance  of  the  experiment,  or  else 
the  results  using  different  limits  will  be  presented. 
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EXTENSION  TO  COMPLEX  SYSTEMS 

Unfortunately,  the  case  of  a  system  consisting  solely  of  non¬ 
identical  parts  in  series  is  relatively  rare.  More  often  one  is 
interested  in  systems  having  a  mixture  of  series  and  parallel-connected 
parts.  The  following  diagrams  illustrate  some  of  the  possibilities: 


The  squares  containing  numbers  denote  parts  of  the  ith  type, 
and  natural  subdivisions  within  systems  may  be  called  components.  For 
example,  the  first  figure  represents  a  system  of  three  parts  in  series 
and  one  component  of  two  identical  parts  connected  in  parallel. 

In  the  following  examples  the  methods  used  are  those  suggested 
for  use  in  more  general  models. 

In  the  first  illustration,  suppose  that  for  the  i^1  part,  n  trials 

have  been  conducted.  Then  from  the  section  on  parametric  bounds, 

3  o 

P  =  TT(1'0(1’(0>  £=2n.q.,  and 

i-1  i=lxl 

(l-qi+2)(l-|/n)  <  P  <  (l-!/3n)5(l-q42). 

Assuming  that  the  numbers  X^  are  available,  it  is  easy  to  find 
two-sided  a-levei  bounds  for  5 ,  and  3-level  bounds  for  q^  based  on 
X^  +  Xg  +  X,  and  respectively.  These  bounds  are  such  that  the 
events,  say  c^  <  I  <  c0>  <  q^  <  d£  independently  occur  with  proba¬ 

bilities  1-a  and  1-3  respectively.  But  if  both  occur  at  the  same  time, 


and  hence 

(l-d22)(l-c£/n)  <  P  <  (l-c1/3n)3(l-d12) 

with  probability  (1-aXl-P).  If  P  is  chosen  small  enough,  then 

(l-d22)(l-c2/n),  (l-c1/3n)5(l-d12) 

-ill  form  a  confidence  interval  for  P  with  coefficient  very  nearly 
1-a.  As  0  approaches  zero,  the  interval  d. ,  d  will  of  course  widen. 

r-  g  Me 

But  i-d^,  l-dg  will  change  by  a  relatively  small  amount.  Hence  one 
may  cnoose  p  very  near  to  zero,  and  the  resulting  confidence  interval 
for  P  will  not  be  appreciably  widened. 

It  may  be  noted  that  one  may  also  lower  a  a  bit  in  order  not  to 
make  the  bounds  so  wide,  at  the  expense  of  a  slight  widening  of  the 
entire  interval,  to  ensure  the  desired  confidence  coeffic..ent.  One  is 
really  free  to  choose  a  and  p  .Just  so  long  as  (l-aXl-P)  attains  the 
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desired  number.  In  this  example.,  the  width  will  be  found  to  suffer 
least  if  a  is  lowered  as  little  as  possible,  and  &  is  made  nearly  zero. 
If  (l-c'Xl-P)  is  slightly  below  the  desired  coefficient,  the  proba¬ 
bility  of  coverage  will  quite  likely  be  higher,  due  to  the  raising  of 
the  actual  confidence  coefficient  noted  previously. 

For  the  second  illustration, 

p  =  TT'^iHi-^^a-q!2)2- 

i*3  1 

1* 

Since  will  generally  be  very  small,  one  may  be  willing  to  ignore  it 
altogether  and  consider 

P  =  Hd-q^d-q^)2, 
i=3 

in  which  case  the  parametric  bounds  are 


(l-q.2)2(l-S/min  n,)  <  P  <  (1-J/max  n  )(l-q2)2 
1  i-3,4  1  i-3,4  1  1 

(assuming  observations  for  the  1th  part  have  been  taken,  and  say, 

5  <  1’.  n,  -  2min  n,  and  min  n4 ,  where  i  *  n,q,  +  n,  q.  ). 
i-3  1  i-3,4  1  i-3,44-  p  p 


If  c1 ,  cg  and  d^  d^  are  a  and  3 -level  two-sided  confidence  limits 
for  5  and  q^  respectively,  then 

(l-d22)c(l-c2/min  ni).  (l-c^max  n1)(l-d12)2 


are  two-sided  confidence  limits  for  P  with  coefficient  greater  than 
or  equal  to  (1-aXl-?)*  If  one  wishes  to  take  account  of  qg,  then 
7-level  limits  b  ,  bg  for  q^  lead  to 

(l-bgi+)(l-d22)2(l-c2/min  n^,  (l-c^rain  ni)(l-d12)2(i-b11 ) 

us  the  limits  for  P,  with  confidence  coefficient  at  least  (l-c:Xl-PXl-7) , 
where  y  is  chosen  to  be  very  nearly  zero. 

In  the  third  illustration, 

P  =  (l-q1)^(l-q^)(l-  [1-  (l-qg)(l-q5) f  ). 

In  this  case  there  are  several  directions  in  which  one  can  proceed, 


some  of  which  are  listed  below: 

A.  One  can  use  the  rather  primitive  inequality 

(l-ql)2(l-q2:i)(l.q3)(l-q4)  <  P  <  (l-q^U  V^l-q^Kl-q^) 

r.ud  find  I'^ur  independent  eonfidenie  limits  at  low  enough  levels  so 
that  the  product  of  the  four  confidence  coefficients  equals  1-a.  The 
limits  will  vary  considerably  in  width  according  to  the  size  of  the 
X^,  and  the  levels  selected.  If  the  q^  happen  to  be  very  close 
to  zero,  this  may  provide  narrow  limits. 

B.  If  for  example 

1-i/min  n,  <  (l-qn) (l-q*)  <  1-S/max  n, ,  that  is 

1=2,3  p  1=2,3  1 

5 

0  <  i  <  rain  (min  n. ,  £  n.  -  2min  n. ),  where  I  »  n„q„  +  n_q  , 

1=2,3  i=2  i*2,3 

then  it  follows  that 

(l-q1),'(l-q4)(i-[|/min  nj2)  <  P  <  (l-U/max  rx±  )2)(l-q1)2(l-qlt). 
i=2,3  1*2,3 

Using  three  independent  confidence  limits,  one  will  in  general  be 
able  to  arrive  at  narrower  confidence  bounds  for  P  than  by  using 
procedure  A. 

C.  If  one  is  willing  to  ignore  terms  in  the  expansion  of  P  which 
involve  products  of  three  or  more  tl  ,  then  one  may  write 

P  =  1  -  (2qi+q1|)  +  qx2  +  2^  -  -  2q2q,  -  q?2 

and  hence 

1  -  (2qx+q4)  -  (q^+q^r  <  P  <  1  -  (2qx+q4)  +  q^  +  Pq^. 

With  arbitrary  one  will  have  difficulty  finding  confidence  limits 
for  (Pq^+q^)  and  (q^+q^).  However  in  the  special  case  that  nx  =  2n4 
and  n.,  =  n2 ,  the  problem  can  be  given  a  simple  solution.  For  in  this 

C.  y 

case  one  may  refer  X^  +  X^  and  Xg  +  X^  to  the  binomial  or  Poisson 

confidence  intervals,  using  the  following  argument:  X  is  distributed 

i  2 

binomially  with  parameters  n^  and  q^  and  variance  n^q^  -  n^q^  .  If 
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one  refers  X  to  the  binomial  distribution  with  parameters  n  /2  and 
*■  o  p 

2q^,  which  has  variance  n^q^  -  2n.(q1  <  n^q^  -  n^q  one  will  be 

making  a  conservative  approximation  in  the  sense  that  (as  in  the  case 

of  the  Poisson  and  binomial  approximations  to  '  the  tail  probabilities 

will  really  be  smaller  than  under  the  approximation.  Thus  a  confidence 

interval  for  2q^  based  on  the  approximation  will  be  expected  to  be 

too  wide.  X^  +  X^  may  be  referred  to  as  binomial  (parameters  n^/2  +  n g, 

l/tn^/2  +  n2]  )  or  Poisson  (parameter  £  *=  (nj/2)(2q^)  +  ngijg)  for  a 

confidence  interval  l’or  § .  But  =  2r.^  means  that  from  any  confidence 

interval  for  |  one  also  derives  one  for  2q^  +  q^  since  5  *  n^Sq^  +  q^) 

and  one  may  divide  the  end  points  through  by  n^.  Similarly  for  qg  +  q^, 

one  may  derive  any  level  confidence  interval  desired.  If  c^,  Cg  are 

two-sided  a-level  bounds  for  2q1  +  q^,  and  d^,  dg  are  two-sided  0 -level 

bounds  for  q„  +  a, .  then  tne  event 
*  'X 

C1  -  2ql+q4  5  c2  and  dl  -  q2+qJ  -  d2 

implies 

'(q2+q5)2  >  -<i22.  q12+2qj*q1+  <  (2q1qJ+)2/2  <  Cg2/2,  and  ( 

l-Cg  <  l-(2q1+qi+)  <.  l-c1> 

and  thus 

i-Cg-dg2  <  P  <  l-c^+cp2/2 

occurs  with  probability  greater  than  or  equal  to  (l-aXl-8) •  This 
procedure  also  will  generally  yield  shorter  intervals  than  procedure  A. 

As  an  illustration  of  how  B(n1>q1)  compares  with  B{n^/2,  2q^), 
the  following  frequencies  have  been  tabled: 


X 

B(100,  .02) 

B(50,  .04) 

X 

B(100,  .02) 

B(50,  .04) 

0 

.13262 

.12989 

b 

.01142 

.01080 

1 

.27065 

.27059 

7 

.00313 

.00283 

2 

.27342 

.27623 

8 

.00074 

.00065 

3 

.18227 

.18416 

9 

.00016 

.00013 

4 

.09021 

.O9OI6 

10 

.00002 

.00002 

9 

.03935 

.03436 

11 

.00001 

.00000 
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0.  Using  the  condition 


4 

2  n.q.  =  e, 
i=l  1 

one  can  find  the  minimum  and  maximum  values  of  P.  These  will  he  func¬ 
tions  uf  5,  and  using  X  one  can  find  confidence  limits  for  these  func¬ 
tions,  and  hence  for  P, 

Unfortunately,  there  are  several  reasons  why  this  procedure  is 
unsatisfactory:  First,  the  bounds  are  not  simple  to  derive  in  many 
cases;  and  second,  their  form  depends  upon  the  particular  range  of 
values  within  which  5  happens  to  lie,  and  in  instances  this  may  be 
quite  uncertain;  and  third,  the  resultant  parametric  interval  is  wider 
than  under  alternative  procedures. 

In  general,  one  can  often  find  bounds  for  P  which  are  almost  as 
narrow,  and  in  many  cases  narrower,  by  considering  groups  of  functions 
of  the  q^  such  as  illustrated  in  procedure  B  above,  which  are  easier 
to  derive.  Terras  of  the  form  TTd-  q^)  are  often  the  main  contributors 
to  the  value  of  P,  while  the  remaining  terms  may  be  roughly  bounded 
with  little  cost  in  terms  of  width  of  the  confidence  interval.  As  to 
which  of  two  or  more  groups  of  parameters  deserve  the  smallest  confi¬ 
dence  level,  trial  and  error  calculations  may  best  provide  the  answer. 

One  can  see  that  with  computation,  values  of  rm  may  be  chosen  to 
minimize  the  length  of  the  interval.  If  one  also  knows  the  approxi¬ 
mate  size  of  the  q^ ,  a  better  choice  of  testing,  procedure  may  be  made. 

In  the  fourth  illustration, 

v  -  TTU-OU-qoOU-i ,4). 

i=l,4  P  J 

Usually  q,  will  be  so  small  as  to  be  entirely  negligible,  and  P  may 

be  taken  to  be 

XJCl-q.  )(l-qpq-.). 
i=l,U  1  ; 

If  one  tries  to  minimize  and  maximize  the  entire  expression, 

subject  to  the  restraint,  the  result,  even  when  all  the  are  equal 
to  rT,  i:.  in  many  sit,w,t  ions  nnr.at  I'-fnotory .  The  bounds  for  one  range  of 


values  ol'  4  =  n^q^  +  are  4/ii,  1,  which,  if  n  is  Large  relative 
to  4,  is  much  too  wide  to  be  of  any  practical  value.  In  the  general 
case,  however,  one  can  often  use  the  bounds  for  (l-q^fl-q^)  and 

to  provide  bounds  for  P.  Pnr  example,  if  f  =  u^q^  +  n^q^, 

T  =•  +  n3q3’ 

(l-f./min  n  )(l-T„2/4n«n,)  <  P  <  (l-Ti/max  n. ) 

1  i=i,4  1  2  2  5  "  "  1  1-1,4  1 

is  valid  when 


T\<  min  (  £  n , -2min  n  ,  min  n. ),  and  T0  <  min  (n0,  n,). 

1  1-1,4  1  i-l,4  1  i-1,4  1  £  d  5 

gives  bounds  for  and  Xg  +  X^  gives  bounds  for  IV, » 
where  the  product  of  the  two  confidence  coefficients  is  l-a.  In  this 
case  a  narrower  interval  will  result  if  the  confidence  coefficient 
for  is  taken  much  closer  to  l-a  than  the  coefficient  for  ?2* 

In  summary,  to  find  a  confidence  interval  for  the  reliability  of 
a  given  system,  the  following  approaches  are  suggested.  Some  numerical 
trial  and  error  may  be  necessary  to  select  the  most  promising  parametric 
interval. 

(1)  The  upper  and  lower  bounds  for  P  are  found  under  the  restric¬ 
tion  £  np4p  =  6.  The  upper  bound  for  one  range  of  values  of  4  will 

be  the  solutions  of  the  equations  dp/dq^  -  Xn^  (i=l, ...,k),  X  /  0,  and 
4  =  £  n^q^j .  The  bounds  for  other  ranges  of  4  may  be  found  from  trial 
and  error  of  various  numerical  quantities  in  P.  The  lower  and  upper 
bounds  will  both  be  functions  of  4  and  the  known  constants.  A  confi¬ 
dence  interval  for  4  based  on  £  X^  can  by  appropriate  algebraic 
manipulation  determine  the  confidence  interval  for  P. 

(2)  Alternatively,  P  is  separated  into  products  of  simpler 
functions,  such  as 

t  TTu-qJ]  [  TTU-Ol* 

i  1  j  0 

and  the  parametric  bounds  for  each  product  found,  as  in  illustrations 
3  and  4.  Then  by  trial  and  error  the  confidence  coefficients  of  limits 
for  the  bounds  of  each  product  are  determined  so  that  the  probability 


of  covering  P  will  be  at  least  1-a,  and  the  interval  will  have  small, 
if  not  minimum  width.  If  there  exists  prior  information  that  certain 
q, ' s  are  small,  then  one  may  be  willing  to  neglect  powers  of  q^  greater 
than  two,  simplifying  the  computations.  With  a  small  amount  of  trial 
and  error  and  even  a  vn  true  idea  about  the  size  of  the  q^,  one  may 
determine  the  various  confidence  coefficients  in  the  product. 

(3)  If  both  (1)  and  (2)  are  unsatisfactory  approaches,  one  may 

expand  P  in  terms  of  the  which  may  then  be  divided  into  homogeneous 

polynomials  of  ascending  dimension.  The  first  group  will  be  a  linear 

combination  of  the  q^  the  second  will  be  a  quadratic  form,  etc.  Often 

the  may  be  felt  to  be  so  small  that  forms  of  higher  order  can  safely 

be  neglected.  At  any  rate,  the  linear  group  will  contribute  most  to 

the  value  of  P,  and  by  the  device  illustrated  in  method  C  of  the  third 

illustration,  one  may  obtain  bounds  for  this  group,  and  rougher  bounds 

for  the  other  groups,  and  hence  for  P,  In  this  case  one  must  have  at 

least  some  of  the  in  certain  known  ratios  to  each  other.  If  the 

experiment  has  already  been  conducted  and  this  is  not  the  case,  one  may 
*  # 

pick  a  new  set  of  s .  such  that  the  n.  are  properly  related  to  each 
other  and  <  ni  for  all  i.  Then  if  one  randomly  selects  obser¬ 
vations  from  the  n.  Bernoulli  observations  previously  obtained  using  a 
1  * 
table  of  randan  numbers ,  the  new  total  number  of  failures  in  n . ,  say 

X^,  will  be  distributed  binomially  with  parameters  n^  and  q^.  One 
may  then  proceed  as  before,  using  the  approach  of  method  C,  illustration 
3,  having  thrown  away  some  of  the  available  information,  however.  This 
method  is  a  variation  of  an  approach  to  the  confidence  interval  problem, 
due  to  Dr.  D.  H.  Evans  of  Bell  Telephone  Lahoratorleu,  appearing  in 
Reference  2. 

(4)  If  none  of  the  previous  attempts  succeeds,  one  may  search  for 
weaker  bounds  for  P  and  use  independent  confidence  intervals,  such  as 
described  in  method  A  of  the  third  illustration. 

Prior  investigation  into  the  desirable  sizes  of  the  ni  in  general 
will  often  simplify  computations,  in  addition  to  giving  a  narrower 


confidence  interval.  Ail  the  remarks  and  examples  given  above  also 
apply  to  the  case  in  which  one-sided  intervals  are  desired.  The  lever 
parametric  hound  of  P  is  usually  easy  to  determine  from  a  small 
amount  of  numerical  trial  and  error. 
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NUMERICAL  EXAMPLES 


In  the  following  compulations  the  Poisson  confidence  limits 

were  used  in  preference  to  the  narrower  binomial  intervals  for  the 

following  reasons:  (l)  For  3 maple  sizes  as  large  as  in  the 

examples,  there  is  little  difference  between  the  two;  (2)  The 

binomial  limits  are  not  easy  to  compute,  it  being  necessary  to  use 

approximating  formulas  for  the  incomplete  beta  function  in  some 

cases;  and  (3)  The  Poisson  limits  are  immediately  available. 

In  problems  where  the  sample  size  is,  say,  less  than  50,  it  is 

advisable  to  use  binomial  limits. 

Consider  the  first  illustration  in  the  preceding  section,  and 

suppose  that  n,  =  ng  =  =  n^  =  500,  and  the  numbers  of  failures 

observed  are  «  3»  Xg  *  1,  X,  =  10,  and  =  2,  and  hence  the  sum 

over  all  X.  equals  l6.  If  one  fixes 
1  4 

i  =  E  n. n. 

i-1  1  * 

for  various  values  and  investigates  the  behavior  of  P,  it  may  be  seen 
that  for  5  <  500, 

1-1/500  <  P  <  l-U/500)2,  where  P  *  f[  (l-q^d-q^2). 

i=l 

In  this  case,  the  .05-level  two-sided  Poisson  confidence  bounds  for 
5noi  are  9*598,  25.400  (Crow  and  Gardner).  Accordingly,  the  confi¬ 
dence  bounds  for  P  are  . 9^92,  . 999^ . 

Alternatively,  one  may  make  use  of  the  inequality 

(l-q42)  (1-5/500)  <  P  <  (l-^/l50'3)5(l-q1+2) 


where 


I  =  L  n.q 


i=l 


iV 


The  .05-level  two-sided  confidence  bounds  for  t  are  8.102,  22.945 
(Crow  and  Gardner)  and  the  .0^1-level  two-sided  interval  for  q^ 
(equal-^ails  Poisson,  using  the  yf  distribution  as  tabled  in 
Bionet rik*  'T’a'hipq  fnr  o+ni  i nt.1  ci ans )  ^  (.o .0345).  The  resulting 


.05-level  interval  for  P  is  (.9530,  .9839)*  appreciably  narrower  than 
before. 

A  third  procedure  in  this  illustration  involves  expanding  P 
in  terms  of  the  q^^  and  ignoring  the  non-quadratic  and  linear  parts. 
The  result  is  that 

i  2 

P  -  1  -  I  +  £  q.q.  -  q4  , 
i<J  0 


and  maximizing 


one  arrives  at  the  inequality 

l-!/500-q4?  <  P  <  l-l/SOO+I^^OO)*. 
k 

In  this  procedure  the  one-sided  .0  1-level  upper  bound  for  q^  may 
be  used  with  the  two-sided  .05-level  bounds  for  I .  The  bounds  for  £ 
are  again  8.102,  22.945,  and  the  upper  bound  for  q^  is  .0335  (using 
the  x2  tables  again '  Performing  the  calculations,  the  .05-level 
bounds  for  P  undt  •  the  above  assumptions  are  .9529*  *9845. 

From  this  example,  it  is  again  shown  that  the  method  of  minimizing 
the  entire  expression  P^,,,,^)  subject  to  £  n^q^  ■»  |  does  not 
necessarily  result  in  the  shortest  confidence  intervals.  For  a 
particular  set  of  values  it  will  be  advantageous  to  gain  an  idea 
of  how  wide  the  resulting  parametric  interval  will  be  under  more  than 
one  of  the  approaches  suggested.  There  is  nothing  "illegal"  about 
picking  the  method  which  yields  the  narrowest  parametric  interval  for 
a  particular  problem. 

There  have  been  various  procedures  suggested  for  finding  con¬ 
fidence  intervals  for  P,  two  of  which  have  appeared  in  Reference  2. 

The  first  method  described  in  the  report  depends  on  the  Tchebycheff 
Inequality,  and  for  the  previous  values  of  X^,  ,  yields  two-sided 

.05-level  bounds  of  .9393,  1.0000.  Assuming  that  the  mode  equals 
the  mean,  the  interval  becomes  .9503,  .9940.  The  second  method,  due 
to  D.  H.  Evans,  depends  on  a  randomization  procedure  and  in  the 


in 


above  illustration  necessitates  using  only  250  observations  for  each 
of  the  first  three  components.  The  number  of  failures  Y  out  of  250 
complete  systems  constructed  ranges  between  0  and  16.  For  a  few  of 
the  values  of  Y  (the  expected  value  of  Y  is  approximately  7)  the 
intervals  are  as  follows: 

Interval  Bounds 

Y  (Poisson;  Crow  and  Gardner) 


a  =  .05 

5 

.968 

.9  71 

7 

.9^5 

.987 

10 

.930 

.979 

15 

.905 

.968 

In  comparing  these  results  with  the  previous  Poisson  limits,-  it  is 
+0  be  noted  that  the  randomization  approach  suffers  from  the  disadvan¬ 
tage  of  being  unable  to  utilize  all  the  information  available  in  the 
samples.  On  the  other  hand,  if  =  ng  =  n^  =  500  and  n)+  =  1000,  one 
may  obtain  a  more  meaningful  comparison  of  the  results  of  the  two 
procedures.  In  this  case,  Y  may  vary  from  iti  to  IS,  and  gives  the 
following  table  of  intervals; 

TABLE  XI 

Confidence  Interval 

Y  P(Y)  (Poisson;  Crow  and  Gardner)  Length 


11 

.0000 

12 

.0024 

.9593 

.9866 

.0273 

15 

.0808 

•  9?  (5 

.9806 

.0293 

14 

.9158 

.9541 

.9838 

.0297 

15 

.0009 

.9525 

.9838 

.0313 

Expected  Length  .0297 

The  previous  Poisson  limits,  using  the  bounds 

(l-q^) (1-1/500)  <  P  <  (l-l/l500)3(l-q42) 

give  the  .05 -level  interval  (.9558,  . 9859) ,  with  length  .0501.  While 
the  randomization  procedure  almost  always  yields  a  slightly  shorter 
interval  (with  probability  .9991);  the  Poisson  interval  is  for  any 
given  data  a  fixed  interval.  Only  the  interval  corresponding  to  the 


most  probable  value  of  Y  Is  completely  included  in  the  Poisson  interval 
in  this  example. 

For  more  complicated  models  it  becomes  difficult  to  calculate 
the  probability  density  of  Y  even  when  the  mo3t  favorable  sample  size 
relations  are  attained,  and  hence  to  compare  the  resultant  limits  with 
the  Poisson  limits.  However  the  mean  and  variance  of  Y  may  be  computed 
and  confidence  intervals  for  values  of  Y  reasonably  near  EY  may  be 
compared  to  the  Poisson  limits  with  respect  to  width  and  centering. 

It  appears  that  when  information  must  be  thrown  away  to  use  the  random¬ 
ization  approach,  the  Poisson  interval  based  on  parametric  bounds  may 
be  appreciably  shorter.  In  the  optimum  case  for  the  randomization 
method,  the  interval  seem3  to  be  slightly  shorter  on  the  average  than 
the  corresponding  Poisson  interval,  at  the  price,  however,  of  an 
increased  variability. 

If  in  the  same  illustration,  n^  =  500,  n 2  =  250,  n-  =>  500,  n4  =  500, 
and  Xx  =  X2  =  Xj  =  =  0,  then  for 

To  <  Z  n.  -  3min  n . ,  that  is ,  To  <  500, 
i=l  1  i»l,2,3  1 

one  can  use 

I-Tj/250  <  P  <  l-(Tj/500)2  or 


(l-T2/250)(l-q42)  <  P  <  (l-T2/500)(l-qp, 


where 


It 

T  =  E  - >  n  ,  T  _ 
'l  ,  'i‘i'  "2 


n«0. 


iV 


i=l  x  -  i=l 

If  a  . 025-level  lower  bound  for  P  is  desired,  one  can  use  the  upper 
bound  in  the  Poisson  equal-tails  .05  confidence  limits  for  T^,  T0. 

For  all  =  0,  the  lower  bounds  for  P  turn  out  to  be  (using  the  .0^1 
upper  Poisson  bound  for  q^)  .9852  and  ,9$t7  respectively.  In  this 
system,  the  first  inequality  will  provide  lower  limits  when  X^  =  0, 
and  she  second,  when  XJ(  is  positive.  As  remarked  before,  it  is 
undesirable  to  select  the  confidence  bounds  on  the  basis  of  the  partic¬ 
ular  s:  ancle  point  observed,  and  hence  the  greater  lower  parametric 


■bound,  attained  by  the  second  interval  (for  q,^  <  .7)  seems  to  make 
it  preferable. 
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of  functioning)  of  the  system  Is  constructed.  of  functioning)  of  the  system  is  constructed. 

In  the  series  case,  vith  k  parts,  the  minimum  and  maximum  or  P  (is  a  function  In  the  series  ease,  with  k  parts,  the  minimi  ana  maximum  of  P  (as  a  functioi 

of  3l’s)  Is  found.  In  the  general  case,  that  is,  for  a  system  having  a  mixture  of  q  «*)■  is  found.  In  the  genera],  case,  that  is,  for  a  systwi  having  a  mixture 

c:  series  and  parallel-connected  parts,  parametric  minima  and  maxima  for  P  are  of  ssries  and  parallel-connected  parts,  parametric  minima  and  rax  in  to-  P  arc 

found  using  several  approaches.  found  using  several  approaches. 


